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Abstract This study reports on experimentally observed near-wall reverse flow events in a fully developed 
flat plate boundary layer at zero pressure gradient with Reynolds numbers between Rer = 1000 and Rcr = 
2700. The reverse flow events are captured using high magniflcation particle image velocimetry sequences 
with record lengths varying from 50,000 to 126,000 samples. Time resolved particle image sequences allow 
singular reverse flow events to be followed over several time steps whereas long records of nearly statistically 
independent samples provide a variety of single snapshots at a higher spatial resolution. The probability 
of occurrence lies in the range of 0.01% to 0.1% which matches predictions made with direct numerical 
simulations (DNS). The self-similar size of the reverse flow bubble is about 30-50 wall units in length and 
5 wall units in height which also agrees well to DNS data provided by Lenaers et al. (ETC 13, Journal of 
Physics: Conference Series 318 (2011) 022013). 
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1 Introduction 

The occurence of near wall flow reversal and with it the presence of negative values of the local wall shear 
stress Tw of turbulent boundary layers have been subject of debate over the past decades. Eckelmann [3] 
postulated that near wall reverse flow was not possible and experimentalists have rarely, if at all, observed 
this somewhat counter-intuitive flow phenomen. On the other hand a variety of direct numerical simulations 
(DNS) suggest the opposite. Eor DNS of zero pressure gradient turbulent boundary-layers (ZPG TBL) events 
of negative shear stress have been reported by Spalart and Coleman m and also for a turbulent channel flow 
by Hu et al. [5]. Similar observations have been made by Lenaers et al. [9] using simulations of turbulent 
channel flow up to Rcr = 1000 and for turbulent pipe flow by Khoury et al. [8]. Cardesa et al. [2] also 
confirm the existence of areas of vanishing wall shear stress in DNS of turbulent channel flow at Rer = 934 
and Rcr = 1834 and associate these so-called critical points with large scale structures that extend up to 
800 wall units downstream. 

Common to the observations of the DNS data is that with increasing Reynolds number both the oc¬ 
curence and the magnitude of the negative axial/streamwise velocities increase. Lenaers et al. [9] report 
reverse flow occurence of 0.01% for Rcr = 180 increasing to 0.06% for Rcr = 1000. In their DNS of fully 
turbulent channel flow Hu et al. [5] report a probability of negative wall shear {jw < 0) of 0.003% at 
Rer = 90 increasing to 0.085% at Rer = 1440. 

Due to their predicted low occurrence reverse flow phenomena have only been observed rather seldom in 
experiments involving ZPG wall bounded flows. To properly catch these events long records are necessary 
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Table 1 Global parameters of the boundary layer experiments with estimated values given in parenthesis 


measurement location 

Xo 

[m] 

3.21 

3.21 

6.76 

6.76 

free stream velocity 

Uoo 

[m s-l] 

5.0 

9.0 

5.0 

9.0 

boundary layer thickness 

s 

[mm] 

(76) 

(68) 

(140) 

(124) 

wall shear rate 

7 = du/dy\o 

[s]-i 

2950 

8300 

2650 

7420 

friction velocity 

11 j- 

[m s-l] 

0.212 

0.355 

0.201 

0.336 

Reynolds number 

Rex 

XlO® 

1.06 

1.90 

2.24 

4.03 

friction Reynolds number 

R^e^f 


1066 

1590 

1842 

2740 

wall unit 


/rm 

71.8 

42.8 

75.7 

45.3 


which until recently has only been possible for single point techniques, for instance through the use of 
laser Doppler velocimetry in a ZPG TBL as reported by Johansson [6]. At the same time the employed 
measurement technique needs to provide adequate spatial resolution as the reverse flow structures observed 
in DNS data are both short-lived and restricted to the viscous sublayer (0(5y~^)). Using the micro pillar 
shear stress imaging technique, Briicker [1] has recently been able to visualize the areas of reverse flow on 
a fiat plate turbulent boundary layer at Rcr ~ 940. 

Flow topology can nowadays be obtained through particle image velocimetry (PIV), yet, in comparison 
to single point techniques, PIV is generally restricted in acquisition frequency, number of samples and 
measurement uncertainty. This can be partially overcome by restricting the camera field of view which 
allows both an increase of sample rate and sample count ca. The following reports on PIV measurements 
in the near wall area of a ZPG TBL using sample counts of up to 120,000 which is shown to be sufficient 
to capture several instance of reverse flow events. 

The PIV measurements were primarily conducted to characterize the upstream conditions for a different 
experiment performed further downstream within the 20 m long test section. Long records, some of which 
are temporally resolved, enable the capture of rare events such as those described in the following. 


2 Wind tunnel facility 

The measurements were performed at the turbulent boundary layer wind tunnel at the Laboratoire de 
Mecanique de Lille (LML). The measurement locations are located about 3.2m and 6.8m downstream of 
the boundary layer trip position which itself is located at the junction between the contraction nozzle and 
the 2x1 m^ rectangular test section. The tripping device consists of a 4mm rod attached to the tunnel wall 
followed by a 100 mm wide strip of coarse sandpaper (roughness ^ 25grid). Full optical access to the 20 m 
long rectangular test section is provided by large glass windows on all for sides. 

Data was acquired at two free stream velocities of Uoo = 5m/s and Uoo = 9m/s with the wind tunnel 
stabilized to within ±0.01 m/s. Temperature stabilization was set at 20.0 ± 0.1 °G. Table[^ provides the 
relevant parameters of the turbulent boundary layer at the specific measurement conditions. While the 
friction velocity can be retrieved directly from the PIV measurements - using the methodology described 
in m - other parameters such as the boundary layer thickness are estimated from theory since only the 
lower portion of boundary layer was captured by the high resolution PIV measurements. 

For the analysis of the data the following orthogonal coordinate system is adopted: the X-axis is aligned 
in streamwise direction, Y is wall-normal and Z is aligned in spanwise direction. 


3 PIV measurement and post-processing 

Following the procedures laid out El only a narrow wall-normal strip was imaged by the PIV camera(s), 
primarily to obtain the wall-normal velocity profile and related higher order statistics. 

Two camera types were used to capture different image sequences. A high-speed GMOS camera with 
36GB of RAM (Dimax-S4, PGO GmbH, Germany) could capture more than 126,000 frames at 6.7 kHz 
to provide continuous time records. By reducing sample rates to 1-2 kHz statistical independence of the 
samples was improved while maintaining a similar sample count. Additional measurements were performed 
using a scientific GMOS PIV camera (Edge 5.5, PGO GmbH, Germany) which featured increased sensitivity 
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Table 2 PIV parameters of the boundary layer experiments 


Camera model 



PCO Dimax-S4 

PCO Edge 5.5 

pixel size 



11.0 X 11.0 

6.5 X 6.5 

magniflcation 

m 

[/xm pixel” 

25.4 

14.1 

image size 

H xW 

[pixel] 

200 X 1008 

200 X 2560 

fleld of view 

w X h 

[mm2] 

5.08 X 25.6 

2.82 X 36.1 

pulse separation 

at 5 m/s 

[fis] 

150 

100 


at 9 m/s 

[^ls] 

100 

65 


and higher spatial resolution. This camera was operated at a double frame rate of 200 Hz to capture long 
records of statistically independent samples. 

The roughly 5 mm wide measurement area was illuminated by a pair of externally modulated continuous 
wave lasers (Kvant Laser, Slovakia) with a combined output power of about 10 W at a wavelength of 520 nm. 
The non-collimated laser beam with a size of about 6x2 mm^ was focussed into a uniform 6 mm wide light 
sheet using a cylindrical lens (// = 200 mm). The resulting waist thickness was on the order of 200/xm 
before entering the wind tunnel glass panel from below. 

Seeding was provided globally in the closed circuit wind tunnel. Consisting of an evaporated-recondensed 
water-glycol mixture, it was introduced in the diffuser downstream of the 20 m long test section just upstream 
of the fan. The size of the aerosol droplets was estimated at l/xm with a lifetime in the order of 10 minutes. 

A telephoto lens (Zeiss Apo-Tessar 300mm/f2.8) with a 100 mm extension tube imaged the near wall 
region with a sufficiently high magnification (m = 0.44) at working distance of ^ 1.1m from the tunnel’s 
centerline. For the high-speed camera with a pixel size of 11 /xm this results in a magnification 25.4/xm in 
object space. For the sCMOS camera, with 6.5/xm pixel pitch, the spatial resolution improves to 14.1 /xm 
per pixel. To make use of nearly the full aperture of the objective lens the optical axis was inclined about 
1.5° with respect to the tunnel wall. 

The acquired data was processed using a conventional 2-C PIV processing package featuring a coarse-to- 
fine resolution pyramid with intermediate image deformation (PIVview2C, PIVTEC GmbH, Germany). To 
obtain reliable mean velocity data and statistics within close proximity to the wall a high aspect ratio image 
sampling window of 64 pixels width and 8 pixels height was chosen. This corresponds to 1.63 x 0.20 mm^ for 
the high-speed camera and 0.90 x 0.11 mm^ for the sGMOS camera. For the latter the sample has an effective 
size of 12.6x^ x 1.57//^ at Uoo = 5m/s increasing to 21. Ix^ x 2.64//^ at Uoo = 9 m/s. The sample overlap 
was set at 75%. Estimates of the mean and unsteady wall shear rate 7 were obtained using a single-line 
cross-correlation approach (i.e. the sampling window only has a wall-normal size of one pixel). 


4 Data analysis 

To verify that the investigated flow is representative for a ZPG TBL the following will describe some relevant 
statistics retrieved from the processed data sets. The measured data is normalized with inner variables using 
the traditional viscous scaling for velocity uf{= Uijur) and length /^(= liUrji')- In this sense the mean 
velocity profiles for u at position x = 3.2 m for both Reynolds numbers are shown in Eig.[^ left. The 
corresponding variances < u'u' >, < v'v' > and covariances < u' v' >, are provided in Eig.[^ right. Both 
plots also contain reference data from a DNS of a ZPG TBL provided by Schlatter et al. Eor the 

most part, the agreement between experiment and simulation is very good (the lines coincide). Due to the 
limited field of view, the outer region of the boundary is not captured. Discrepancies can be observed very to 
the wall for the high-Re case, stemming from the loss of resolution proportional to the reduction in viscous 
length scales. The covariance < uv' > (lower fine grouping in Eig.[^ right) does exhibit some differences 
in both the buffer layer (around 10y~^) and outer region (> 200y~^) for which there currently is no obvious 
explanation. 

Making use of the high spatial resolution near the wall, both the mean and unsteady wall shear rate 
7 = dujdy and with it the corresponding wall shear stress Tw = ydu/dy can be directly estimated from 
the velocity gradient at the wall. The reader is referred to [13 for details on the processing scheme. Fig.|| 
provides representative probability density functions (PDE) of the wall shear stress Tw for several different 
sequences and show good agreement with data published in literature (e.g | 4 l[ 5 ir 7 ll^[T()] . The vertical fine at 
-2.3 marks the position for Tw = 0. When plotted in log-linear form the PDE of the shear stress exhibits 






4 


Christian E. Willert^ 




Fig. 1 Profiles of the mean streamwise velocity (left) and Reynolds stresses for two different free stream velocities obtained 
X = 3.2 m downstream of the tripping device scaled with inner variables. Symbols only shown for < 10 for clarity. 



-2 0,2 4 6 


Fig. 2 Probability density functions of wall shear rate 7 = 
line at -2.3 represents 7 = 0 



for different data sets obtained at Rcr = 1590. Vertical 


several instances of negative shear stress with a probability of less than 0.01%. Since the wall shear stress 
is directly related to the near-wall velocity, the PDFs of the streamwise velocity u at wall distances of ly^ 
and are provided in Fig.[^ Here it can be observed that the reverse flow seems to only appear very 
close to the wall while it is practically absent outside of the viscous sublayer for y^ <b. 

At this point it should be noted that the data in the tails of the PDFs have an increasing likelihood of 
being affected by measurements errors (outliers) rather than representing reliable measurements. Therefore 
the underlying data sets require separate verification to determine a given datum’s validity which, taking 
into account the rather low probability of less than 0.01%, is feasible through visual inspection. This can 
be achieved through velocity-vs-time plots such as shown in Fig.[^ This image is compiled by extracting a 
single column of data from each PIV data set and placing the columns side-by-side such that the resulting 
image has a width of up to 126,000 pixels, depending on the length of the complete sequence. Therefore 
each line in the image represents the velocity record for a given wall distance. 

Reverse flow events can be easily detected by highlighting negative velocities in images such as Fig.[^ 
and retrieving the corresponding single data sets from the sequence for closer inspection. One such event is 
the white spot near the middle of the bottom edge of Fig.[^ The spot has a size of about 6 pixels height and 
12 pixels width, the former giving an indication on the vertical height of the reverse flow bubble (i.e. about 
300/im given a grid step size of 2 pixel). In the present sequence of 126,000 images only two such reverse 
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Fig. 3 Probability density functions of velocity at wall distances of 0.5y+ and 5y+ for Rcr = 1070 



Fig. 4 Time trace of streamwise velocity u' at t/oo = 5 m/s {Rct = 1840) covering 0.9 s (top) and 0.075 s (500 samples, 
bottom). The white region near the center indicates a single reverse flow event. Vertical axis represents wall distance 
0 < y < 95+ (?^ 7 mm) 


flow events can be detected. Taking into account the duration of about 10-15 samples per event results in 
a probability of about 0.01%. 

As the image sequences are temporally well resolved, the evolution of a specific reverse flow event can be 
observed within the narrow field of view. In this sense the particle tracks, compiled through the summation 
of several images, visualize the shape of the flow structure (Fig.[^. The passage of the reverse flow region 
through the field of view is further visualized by a sequence of velocity fields in Fig.for which only every 
fourth frame is shown. The red patch is this sequence defines the zone of negative streamwise velocity u. A 
magnified view of a single data set from Fig.is provided in Fig.[^ Finally snapshots of several different 
’’separation bubbles” are shown in Fig. |8]for a free stream velocity of Uoo = 9 m/s at measurement location 
X = 3.2 m. 


5 Results and discussion 

In total eight sequences at two measurement locations and two free stream velocities were investigated for 
the presence of reverse flow. While some sequences were affected by low seeding density and limited spatial 
resolution, all sequences showed clear evidence of reverse flow in the form of particles moving upstream for 
a certain duration. The appearance of the flow features is on the order of 0.01%. 

In most cases the reverse flow region has a vertical dimension of and a length of about 30 — 40^::^ 
which corresponds to the DNS provided by Lenaers et al [9] (e.g. see Fig. 9 in their publication). Time- 
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Fig. 5 Multi-exposed images showing a near wall separation event in a ZPG TBL at Rct = 1840 

{x = 6.8 m, Uoo = 5 m/s). Full image view at left with detail on left. The tunnel glass wall is located at y = 0. 

Vx [m/s] 


0 12 3 4 



Fig. 6 Sequence showing the evolution of a flow reversal event (red area) in a ZPG TBL at {Rct = 1840). Red patch near 
the wall indicates reversed streamwise velocity. Temporal separation between frames is 600 ys. PIV sampling window of 
32x8 pixels with 75% overlap, vectors down-sampled 4x horizontally, 2x vertically 


resolved sequences show that the ’’separation bubble” traverses downstream through the field of view at a 
convection speed Uc of about 0.1 Uoo or Udur ~ 2.5 (estimated from Fig.[^. While not captured through 
the present measurements, the corresponding DNS indicate that the structures have a spanwise dimension 
on the order of 30^^, that is, their xz shape is roughly circular. 

Due to their rather limited number in comparison to the number of available samples a more detailed 
statistical analysis is not reasonable. Furthermore, the boundary conditions of the wind tunnel have yet 
to be fully characterized. While the wind tunnel test section generally has a slight favorable pressure 
gradient (personal communication), a model was installed at x > 10 m further downstream of the present 
measurement locations. The presence of the model will undoubtedly have some upstream influence on the 
local pressure field which in turn can affect the flow statistics and with it the probability of flow reversal 
events. At the very least, the present investigation shows that the topology of low-probability reverse flow 
events can indeed be captured using sufficiently long PIV data sets. 


6 Summary and outlook 

Through analysis of long PIV data sequences rare events such as small-scale near wall flow reversal could 
be documented at two measurement locations and two free stream velocities of the ZPG TBL. Both the 









































Near-wall flow reversal in a ZPG turbulent boundary layer 


7 



-20 -10 0 10 20 
+ 

X 

Fig. 7 Detail view of center frame of Fig.|^ shown in viscous scaled units (vectors down-sampled 2x horizontally) 



Fig. 8 Various instances of near wall reverse flow indicated by red patches (Rer = 1590). 


probabilty of occurence as well as the shape of the observed reverse flow structure structures agree with 
previous DNS by Lenaers et al [9]. In the present measurement configuration the spanwise extension of the 
structures could not be measured. Multiple-camera (photogrammetric), time-resolved techniques such as 
tomographic PIV [TT] or 3-D PTV [12], or digital holography [15] are ideal candidates to capture the fully 
resolved velocity field of the small separation bubble. 

The present data sets also exhibit strong wall-normal velocity events very close to the wall but have not 
been studied in detail in the present study. 
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